Nonlinear nonlocal diffusion of magnetic flux in thin type-II superconductors and 

Josephson junction arrays: Exact solutions 
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An exact solution of the nonlinear nonlocal diffusion problem is obtained that describes the 
evolution of the magnetic flux injected into a soft or hard type-II superconductor film or a two- 
dimensional Josephson junction array. (The magnetic field in vortices is assumed to be perpendicular 
to the film; the electric field induced by the vortex motion is proportional to the local magnetic 
induction; flux creep in the hard superconductors under consideration is described by the logarithmic 
U(j) dependence.) Self-similar flux distributions with sharp square-root fronts are found. The fronts 
are shown to expand with power law time-dependence. A sharp peak in the middle of the distribution 
appears in the hard superconductor case. 

PACS numbers: 74.60.Ec, 74.60. Jg 
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There are two kinds of problems of magnetic flux evo- 
lution in the type-II superconductors. In the first type 
problems the flux lines are parallel to the surfaces of 
the superconductor plate. In this case the motion of 
every vortex is determined by it's nearest environment, 
so local nonlinear diffusion takes place. Jt-,can be de- 
scribed by nonlinear differential equationafu, which are 
similar to classics] equations for nonlinear diffusion (see, 
for example ,E2l"L3). In the second type problems the 
flux lines are pemendicular to the surface of a thin 
superconductorli3~E3. Then every vortex is influenced 
by all vortices in the superconductor, and one observes 
a so called nonlocal nonlinear diffusion of the vortices, 
which can described by a complicated integral differen- 
tial equation. Note that the geometry of a fiat sample 
in a perpendicular magnetic field is realized in most ex- 
periments which study the space-time image-pf the mag- 
netic flux in type-II superconductors (seeE2l and refer- 
ences therein). Very similar to that is the problem of the 
magnetic flux penetration into two-dimensional Joseph- 
son junction arraysE3. Macroscopic dynamical equations 
for the magnetic flux in the last case are the same as for 
type-II syperconductor films. 

Up to now the last geometry was theoretically stud- 
ied only in the particular case of the flux flow resistivity 
independent of the magnetic field. The situation where 
the electric field E depends linearly on the current den- 
sity jliSEij was considered numerically. Switching on of 
a constant magnetic field in the situation, in which the 
electric field depends exponentially on j and the flux flow 
resistivity is again independent of the local magnetic in- 
duction B, was treated exactly in paperEH. (Note that, 
strictly speaking, one can introduce the quantity "flux 
flow resistivity" only if E cx j. When we use this term 
for brevity, we mean that E(B,j) can be written in the 
form E(B,j) = p(B)e(j), where p and e are some func- 
tions.) 

In the present paper we shall exactly describe the evo- 
lution of the magnetic flux injected into the film of a 
type-II superconductor, in which the flux flow resistivity 



is proportional to the local magnetic induction (that is 
the most widespread situationEj) . The cases of other field 
dependencies of the flux flow resistivity will be discussed 
only briefly. Note that the flux is assumed to be injected 
into a middle of the film instead of penetration of the 
magnetic flux from boundaries which is studied usually. 

Let us consider the evolution of the magnetic flux in- 
jected in the infinite thin film (the xy plane) of a type-II 
superconductor (the flux lines are perpendicular to the 
surface). We assume the problem to be homogeneous 
along y, so the local magnetic induction B depends only 
on the coordinate x and on time. The current flows along 
y. An applied magnetic field is absent. We write out im- 
mediately the equation for B{x,t) in the most general 
form (afterwards we shall show how it can be obtained): 



\B(x) 



du 



dB(x) =D d 
dt dx 
B{u) 



sign 



da 



B{u) 



(1) 



Here we have introduced the dimensional coefficient D 
and the constants p > 0, q > 0. It is convenient for us to 
use Eq. (FT) rather than the equation for the current den- 
sity p(seel3ll3) or than the equation for the electric field 
(seetll). 

Let us first show how Eq. (]l|) can be obtained for a 
soft type-II superconductor. The continuity equation for 
the vortex density n (B = &on, $o is a flux quantum) 
has the form: 



On 
~dt 



divJi 



0. 



(2) 



In the simplest case, the vortex flow J„ which is propor- 
tional to the electric field, can be expressed in terms of 
the vortex density and the vortex velocity: J v = nv (v 
is the velocity of vortices; the sign of n(x) is plus if the 
magnetic field in vortices in the point under considera- 
tion is parallel to the z axis and it is minus otherwise; 
we assume that vortices and antivortices annihilate in- 
stantaneously) . One can write for a soft superconductor: 
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rjv = (<f>o/c) j, where j is the current density, c is the ve- 
locity of light, and the viscosity r) is related to the normal 
phase resistivity p n by the equation Qq/v — PnC 2 /H C 2 
(H C 2 is the upper critical magnetic field). Thus we have 
for a soft superconductor in the considered geometry: 



dt en dx 



(3) 



This equation is also available in the case of the flux lin^ 
parallel to the surface of a soft superconductor - see 
where we discuss the reasons for the modulus of B to 
appear. However, in this situation j = — (l/AncjdB/dx, 
so 



dB 
~dt 



$o d 
4ircn dx 



B 



dB 

dx 



(4) 



In a thin infinite superconductor film of the thickness d, 
B(x, t) related to the current density j(x, t) by Ampere's 
law, which reads 



c 7_™ u — x 



(5) 



i.e. the Gilbert transformationE3. It may be inverted to 
give 



2n 2 d 



du 



B{u) 



(6) 



Inserting the expression (Q) for j into Eq. one ob- 
tains Eq. CjTJ> with q = p = 1 and D = $ / \2pdr]). Ojhe 
corresponding stationary problem was considered ir£j.) 

Frequently, in the flux creep regime, one can assume 
that the effective creep activation barrier grows loga- 
rithmically with decreasing current: U(j) = t/o In |j c /j|, 
where Uq is the characteristic scale of the pinning en- 
ergy barrier, j c is the critical current density, so the 
vortex velocity v = luUh \ j / j c \ U °^ kT j / jc, where and 
ujh are the averaged hopping distance and the hopping 
frequency correspondingly (sedd and references therein). 
Using again Eqs. (||) and (^), we obtain Eq. (||) with 
q = l,p= l + U Q /kT,D = l h uj h (c/2ir 2 dj c ) p . If before the 
injection of the flux, there is the homogeneous magnetic 
induction B$ in a sample, then for the extra magnetic 
induction SB(x, t) = B(x, t) — Bo <§; Bo one obtains Eq. 
(Q) with q = 0, so real linearization is absolutely impos- 
sible for hard superconductors. 

Let the magnetic flux $ = J dxB(x) > be injected 
in the sample at the initial moment. (We assume that 
B{x) is nonzero in a restricted range around the point 
x = 0.) At long times one may search the solution in the 
following self-similar form: 



B(x,t) = C x t- b 



C 2 t a 



(7) 



where C\ and C2 are constants. Inserting this form into 
Eq. (nh and comparing powers of t in all terms of the 



equation, we get the following relation for the exponents 
a and f3: 



+ P(q+p-l) 



(8) 



Using the condition of the flux conservation $ = 
J dxB{x,t) — const, we get a — /3 = l/(q + p). Thus 
solution to Eq. (|l|) can be written in the following self- 
similar form 



B(x,t) 



tD 
IT 



-V(9+P) 



\($q+p-ltD) 1/{q+P) 

where 6(£) is defined by the equation 
1 



(9) 



dt 



■ p dt, dt, 

KO 



sign 



dc, 



HO 



(10) 



Note that a full derivative appears in the right hand part 
of Eq. (|l(]) because of the equality of the exponents a 
and (3 in the regime of conservation of the injected flux. 
Thus, for the values of £ where 6(£) is nonzero, we get 
the equation 



HO 



i/p 



b {1 - q)/p (0m 1/P sign(0. 

(11) 



(We assume that b(— £) = &(£), so the integration con- 
stant appears to be zero). j— . 

Using known analytical properties of Cauchy integralcHl 
near the edge of a contour and the symmetry of &(£), one 
may see from Eq. ml]) that 



6(0-6(0) = — 



1 



7T \q + p 



l/p 



6 (i-?)/P( ) cot — 
2p 

(12) 



when £ — > (here p > 1/2). If the function 6(£) has tails 
at infinities, then one obtains immediately from Eq. (fTA) 



b(0 = (q + p) 1/(1 - q) 



dCHO 



P/(l-«) 



lei 



-(p+i)/(i-«) 

(13) 



at |£| — > co. 

If q < 1 (including the physically interesting case 
q = 0), there are the feature ( |l2"| ) at the center and the 
tails (O) at |£| — * 00. (Note that d(b(Q < 00 since 
(p + 1)/(1 — q) > 1 for q < 1.) In particular, when 
q = 0,p = 1, (Q) and (^lj) turns to be linear equations, 
and ( |ll| ) has a simple solution 

const 



HO 



(14) 
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The relation (|l3|) is obviously impossible if q > 1, and 
in this case &(£) may be nonzero only near the point 
£ = 0, so 6(0 = if |£| > &> > I- When q > 1, we 
shall search the solution among such functions. It is con- 
venient to use the function h(£) which is nonzero in the 
range — 1 < £ < 1, instead of 



K0 = Co 



i/(g+P-i) 



(15) 



The parameter £o will be found at the end of our calcula- 
tions. The integral in Eq. ( pd| ) becomes the integral be- 
tween the limits — 1 and 1 . Using known analytical prop- 
erties of the Cauchy integral near the edges of the contour 
of integration, we see from Eq.(|i"T|) that h(£) oc y/l — £ 2 
near the points £ = ±1 when q = 1. h(£ ~ ±1) can also 
be found in the case q > 1. In particular, 



M0 = 



i 



l/(p+9-l) 



x In 



g-1 



p/(p+«-i) 



1-C 



-p/(p+<?-i) 



(16) 




FIG. 1. Scaling functions /i(£) determining the evolution 
of the magnetic flux injected into a type-II superconductor 
film (see Eqs. (§), @, @, and ([§). Curves 1,2,3 cor- 
respond to: g = p = 1 - the case of a soft superconductor; 
q = l,p = 3; and < g < oo,p — > oo - the case of a hard 
superconductor. 

Figs. [j] and || show scaling solutions h(£) for g = 1 and 
different p. The case p > 1 (Fig. [l]) is interesting for us. 
The solutions with two maxima for < p < 1 (Fig. 
^J) are shown here for generality. In particular, in the 
case of a soft superconductor {q = p = 1) 



when 1 — £ <C 1. We do not discuss this case in detail. 

Let us consider the most interesting case q =p-il, for 
which it appears possible to find an exact solutionis. One 
may invert the Gilbert transformation, using relations 
for functions restricted at the edges of a finite intervalEa. 
Then we obtain 



H0 = 

When q — 1 , p — > oo 



2tt 



K0 = 



In 



Co = tt/2, 



(19) 



(20) 



MO 



p+ 1 



1/p 



i i 



|C| Vp ^n(C) 



TT 2 p \p + 1 



1 3 >9 

— 1- -• 1 — £ 2 

2p' '2' ? 



i/p 



(17) 



where £?(, ) is the beta-function and 2-F1 is the hyperge- 
ometric function. h(0) = ^{p + 1)-Vpb U, . Now 
we recall that $ = j dxB(x, t) — const. Thus, 



£ 0=7r P/(p+i) (p + 



1)B -p/(p+D fl i) 



(18) 




FIG. 2. Same as Fig. |l| but curves 1,2,3 correspond to: 
q — p = 1; q — l,p = 0.6; and q = l,p = 0.11. 



The relations @, (15), (p"7|), and (18) give an exact solu- 
tion of the problem under consideration in the case q = 1. 



and we obtain the solution for the case of the hard 
superconductor at low temperatures. This limit has the 
same form also for other values < q < 00. Here, in each 
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point of the distribution, the integral J d£h(Q /(£ — = 
const, so we are concerned with a nonlocal analog of the 
sandpile problem. As we have noted above, at q = 1, the 
solution has a square- root form near the fronts. Now we 
may see from Eq. ([l7]) that 



h(£~±l) = —(p + l) 
ir z p 



-i/p 



B 



1 1 
2' ~2~p 



v/r^e (2i) 



Now let us show how the current density distribution 
evolves. When q > 1, one may find from the relations 
(PD, (p|), and (pi) that 



tD 



2n 2 d V $ 



i/(g+p) 



+p-:ur)\ 1 /(9+p) 



(22) 



where g(£) = -/VCMOAC ~ 0.; When |£| < 1 and 

9 = 1, .9(0 = OP + l)- 1/p |e| Vp «5n(0- If 9 - 1 and 
> 1, using (0) one gets 



"«>=> +i >- < * +i, " s (H 

1/11 _L_ 3 1_ 1_ 

x e v 2 ' 2 + 2 p ; 2 + 2 p ; ^ 2 



(23) 




FIG. 3. Scaling functions g(£) (see Eq. (|23|) for the current 
density distributions corresponding to the flux distributions 
shown in Fig. [I]. Curves 1,2,3 correspond to: q = p = 1; 
q — l,p — 3; and < q < oo,p — > oo. 

(see Fig. |). Near the points £ = ±1 (at |£| > 1), 

5(C)-5(IC| = 1) = 

"^(P + !)" 1/PS Q> ^) N/iiFl«fln(0. ( 2 4) 

In particular, if g = p = 1, then g(£) = £/2 for |£| < 1 
and 



the limit p —* oo,q < oo, <?(£) = sign{£/) for |£| < 1 
and 5(£) = (2/7r)arcsin(l/0 for |£| > 1. In the lin- 
ear case q = 0,p = L, as can be seen from Eq. (H), 
the current density (|6J) is determined by the integral 
-L OO oo ^(0/(C-0 = ^n S ^/(£ 2 + 7r 2 ). 

Remarkably, in the studied problem, nonspecified ini- 
tial distributions of the magnetic flux evolve to unified 
space-time distributions @ at long times. It is signifi- 
cant here that $ ^ 0. If total flux of the initial distribu- 
tion is zero, then the distribution evolves to a different 
scaling form at long times. 

The duration of the transition from the initial distribu- 
tion of the injected magnetic flux to the scaling solutions 
(see Fig. |l| and Fig. |^) depends on the specific form of 
the initial distribution. Obviously, in the case q > 1, the 
initial distributions should be localized to obtain scaling 
solutions at long times (see the discussion after Eq. Jl^). 

In conclusion, let us compare the obtained solutions 
with the solutions of the local equation 



9B =D d_ 
dt dx 



|./J|'' \<)B/(>.r\ r sign ( ^ 



(25) 



describing evolution of the magnetic flux injected into a 
type-II superconductor plate when the flux lines are par- 
allel to the surfaces. (We assume that the problem is 
homogeneous along a surface.) Let the flux be injected 
insite of a plate. In this case, until the vortices reach 
the boundaries, B(x,t) = Cit-Vta+arift (x/C 2 t 1 ' {q+2 ^) 
where C\ and C 2 are constants, and one can see that 

HO oc [l-(e/0) (1+p)/p ] P/(9+P_1) . This function has 
no such striking peculiarities at the center as in Figs. |l| 
and |. 

We consider above the homogeneous along y axis prob- 
lem. One can also obtain the exponents a and /3 intro- 
duced in Eq. (Q) for the situation, in which the initial 
distribution of the injected flux is localized not along a 
line but in the vicinity of some point on a film, so instead 
of the coordinate x now one can introduce the radius r. 
Then, as it can be shown, in the case of nonlinear non- 
local diffusion in a film, a = (3/2 — l/(2q + 2p — 1) (it 
should be q + p > 1/2). For a parallel geometry (i.e. for 
nonlinear local diffusion in a plate) one can also search for 
scaling solutions with axial symmetry. In this problem, 
a = 0/2= l/(2q + 3p- 1), 2q + 3p> 1. 

As it can be easy understood, the problem described by 
Eq. (Q) is a direct generalization of the evolution problem 
for charge distribution placed on a conducting plane with 
dissipative transport. Analogical problems for charges in- 
jected into a wire and for evolution of yz homogeneous 
charge distributions in 3D can be also formulated. The 
equation for the first one is of the form 



dn (x) d 



du 



n (u) 



\u — x\{u — x) 



(26) 



2 £ fl + \/l £ 2 ) f° r 10 — 1- I n and the last one is described by the following equation, 
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9n (x) d . 



dun (u) sign(u — x) > , (27) 



where n(x) is the charge distribution. The scaling so- 
lutions of Eqs (^6|) and ( |27j ) are of the form @ with 
a = (3 = 1/2. The scaling function &(£) looks like 



(28) 



for the solution of Eq. (Pq), where is a step function, 
and 



(29) 



for the solution of Eq. ( p7f ) . One may compare this func- 
tions with answers for Eq. (|l|) shown in Figs [j] and [|. 
Note that an exact solution can be found generally for 
the equation of the form: dB/dt = d{\B\ F{B}) /dx, if 
the condition F{cB} — c p F{B} is satisfied (c = const ^ 
0,p = const > 0). 

In summary, we have considered problem of nonlin- 
ear nonlocal diffusion of the magnetic flux injected in 
an infinite thin type-II superconductor. We have solved 
it exactly in the most interesting case of flow resistivity 
proportional to the local magnetic induction B. The ob- 
tained flux space-time distributions are of the self-similar 
form with rather striking scaling functions. Two ques- 
tions remain open: how can the obtained distributions 
be observed experimentally, and how do the edges of a 
thin strip affect our solutions? 
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